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Abstract 

A discussion of discrete Wigner functions in phase space related to mutually unbiased bases 
is presented. This approach requires mathematical assumptions which limits it to systems with 
density matrices defined on complex Hilbert spaces of dimension where p is a prime number. 
With this limitation it is possible to define a phase space and Wigner functions in close analogy to 
the continuous case. That is, we use a phase space that is a direct sum of n two-dimensional vector 
spaces each containing p'^ points. This is in contrast to the more usual choice of a two-dimensional 
phase space containing p^" points. A useful aspect of this approach is that we can relate complete 
separability of density matrices and their Wigner functions in a natural way. We discuss this in 
detail for bipartite systems and present the generalization to arbitrary numbers of subsystems when 
p is odd. Special attention is required for two qubits (p = 2) and our technique fails to establish the 
separability property for more than two qubits. Finally we give a brief discussion of Hamiltonian 
dynamics in the language developed in the paper. 
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I. INTRODUCTION 

In a study of thermal equilibrium of quantum systems j26[ , Wigner introduced the famous 
function that now bears his name. There is an extensive literature on the Wimer function 
fo.. co.t„u,o„s va.ab.es The Ute^atu. o. discrete W.gne. function, is .ess exte.s,ve. 

but the importance of discrete phase space in quantum information has revived interest in 
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2J]. In particular, the paper by Gibbons, et. al. contains a useful list of 



the subject 
references. 

In this paper we present a discussion of discrete Wigner functions in phase spaces related 
to mutually unbiased bases (MUB). Our approach differs from the geometric method of 
Wootters in being more operational and closer to the methodology of the continuous case 
^, , but our approach also requires mathematical assumptions which limits it to systems 
with density matrices defined on complex Hilbert spaces of dimension where p is a prime 
number. With this limitation it is possible to define phase space and Wigner functions which 



mimic the continuous case, 
dimensions, see for example 



here does not seem to be any simple way to do this for other 



. A useful aspect of this approach is that we can relate 
the separability of density matrices and their Wigner functions. We discuss this in detail 
for bipartite systems and present the generalization to arbitrary numbers of subsystems. As 
an application of our analysis, we show that for p an odd prime, with a particular choice 
of "phase" parameters, Hermitian operators used in 8] for n p— level systems are tensor 
products of opeators for the individual p— level subsystems. 

The paper is organized as follows. We first briefly review the definition and properties of 
the Wigner function for continuous variables and list the most important properties that are 
retained in the discrete case. Our discussion of the discrete Wigner function makes extensive 
use of generalized spin matrices which are defined in the section III for a singl particle. In 
order to determine a suitable choice of phase space, we are led to consider mutually unbiased 
bases, and this is done in sections IIVI and IVH and further discussed in Appendix IX Bl The 
discrete Wigner function for a single particle is then defined and its properties discussed in 
section IVl The generalization of our discussion to more than one particle begins with section 
IVII The transition to the general case is aided by using the geometry of discrete phase space, 
which is summarized in Appendix IX El In section IVII Dl we generalize the Wigner function 
to dimension p^, and in IVIIII to p". 
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The problem of separability when p = 2 requires special treatment, and in section IVTTl the 
case of two qubits is analyzed. The generalization to more than two qubits appears to be 
impossible by the present technique, this is discussed in section IVlllI Finally, in section IX 
a brief discussion of Hamiltonian dymanics is presented and a simple example using MUB 
is given. Various background and technical issues are discussed in the appendices, including 
the positivity of the density matrix. 



II. WIGNER FUNCTION FOR A PARTICLE MOVING IN ONE DIMENSION 

Let p be the density matrix for a particle moving in one dimension, and let Q and P 
be the position and momentum operators for the particle. We set h = 1 so the Heisenberg 
commutation relation is [Q, -P] = H- It is convenient to introduce the Wigner function as 
the Fourier transform of its characteristic function x, defined by 

Xp{u,v) =tr[pD{u,v)] (1) 

where D is the unitary translation operator 

D{u, v) = e-^("-P-^«) = e-^"^e™^e^™/^ (2) 

These operators form a projective group called the Heisenberg- Weyl group 25j. It is easy 
to show that 

D{u, v)D{a, b)D{u, v)^ = e'^"'^^"^'"'''^ D{a, b), (3) 
where the phase factor is the symplectic product of the operator "indices" , 



(a, b) o [u, v) = bu — av. (4) 



The Wigner function is defined by 



OO POO 
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-OO J — OO 

OO POO 



du / dvtr[pDiu,v)]e-'^''''-P''\ (5) 



To see that this agrees with the standard definition let us compute the trace in the last 
equation using a complete set of eigenvectors of Q, 

Wp{q,p) = J duj dv j dx(s|p|a; + M)e™^e*™/2e-i(9''-P«) 



where Eq. Q was used with 

e~'''^\x) = \x + u). 

Doing the v and x integrals gives 

WpilyP) = — f du I dx{x\p\x + u)6{x -\ q)e 

271 J J 2 



1 

2^ 



du{q-u/2\p\q + u/2)e'P'^. (6) 



The definition of the operators D{u,v) is not unique. There is some freedom in the choice 
of phase, referred to as gauge 

is the standard one, the issue is not so simple for the discrete case. 

Many of the standard properties of the Wigner function can be deduced readily from 
Eq. ®: 

1. the mapping p — > Wp is convex linear. 



2. Wp is normalized, i.e. 



dq J dpWp{q,p) = 1 



which follows from Xp(0, 0) = trp = 1, 
3. W is real since x*p{u,v) = Xp{—u, —v) 
4 . if p' = D{a, bypD{a, h) then 



^pi{q,p) = J^f^ J J dvtr[pD{a,b)D{u,v)D{a,by]e~ 
= j du j dvtr[pD{u,v)e'^''^-''^^]e 



^i(ub—va)~\ ^—i{qv—pu) 

= Wp{q + a,p + b), 
5. the marginal distributions are probability densities, 

dqWp{q,p) = {p\p\p) 
dpWp{q,p) = {q\p\q). 

More generally, if we integrate along a line in phase space we get a probability density 
dq / dpWp{q,p)5{qcose + psinO - qo) = {qo]9\p\qo;e), 

-OO J —OO 

where |go; ^) is the eigenvector of Qe = Qcos9 + Psin6' with eigenvalue go- 



■OO 
OO 
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Finally, to show that the Wigner function is equivalent to the density matrix, we write 
the density matrix in terms of the Wigner function. This is done easily by taking the inverse 
Fourier transform of Eq. (jUj) 

It follows from this equation that 

/QO POO 
/ dpdqWp,{q,p)Wp^{q,p), 
■oo J ~oo 

which is just Plancheral's theorem. 

Proving that a given function W{q, p) corresponds to a density matrix comes down proving 
that the inverse formula leads to a p which is positive (cf. ref |l6|). 

Finally we note that we can define a Wigner function, Wa-, for any operator A for which 
Eq. (0) is defined. 



III. GENERALIZED SPIN MATRICES 



We briefly review some facts about the generalized spin matrices which are of interest 
here and introduce some notation that will be used throughout the paper. We shall use 
letters j, fc, s, t to denote elements of = {0, 1, ■ ■ ■ , d — 1}, the integers modulo d. Let Hd 
be a d -dimensional complex Hilbert space, and let {\k)^ k G Zd} be an orthonormal basis 
of Hd- Let Md be the vector space of complex d x d matrices that act on Hd- This space is 
a d"^ -dimensional Hilbert space with respect to the Frobenius or trace inner product 

{A, B) = tr [A^B) (7) 

for A,Be Md- The set of matrices {|j)(A;|, j,k G Zd} is an orthonormal basis of Md- Let 
1] = r]d = e*^'^/'^, and define the generalized spin matrices as the set of unitary matrices 

d-l 

S*^.^ = ^r/^"^|m)(m + A;| (8) 

m=0 

where index addition is to be understood to be modulo d. This set of d"^ matrices, including 
the identity matrix / = 5^0, forms an orthogonal basis of Md jiol . 
It is not difficult to show that 

Sl = V^'S^,,., (9) 

Sj,kSs,t = V^"^ Sj+s,k+t- (10) 



From Eq. (fTUjl it follows that Sj^k and Sg^t commute if and only if the symplectic product 
(j, k) o (s, t) = 0, where 

(j, k) o (s, t) = ks — jt mod d (11) 
which should be compared with Eq. (@]). We also will need the relation 

cm _ m{m-l)jk/2q /-|r,\ 

The spin matrices can be generated from two matrices: 6*1^0 which is diagonal, and 6*0,1 
which is real and translates each state to the next lowest one modulo d. One can check 
that Sj,k = S{ qSqi- These spin matrices can be viewed as translation operators in a manner 
similiar to the D{u,v) operators for the single particle discussed in section |n] The analog 
to property 4 is 

Since the matrices {'^Sj,k} form an orthonormal basis on the -dimension Hilbert space 
Mrf, they satisfy the completeness relation 

^ d-l 

- ^ S,,k tr{SlA) = A, (14) 

j,k=0 

where A G M^. This set of spin matrices has appeared repeatedly in the mathematics and 
physics literature, for example j^ls, 11, 19, 21, 27 1 among others, and is often also referred 



to as the (discrete) Heisenberg-Weyl group. 

Finally we define a set of orthogonal one-dimensional projection operators that we will 
need. Let p be a prime number. For (j, k) 7^ (0, 0) and < r < p — 1 

where ^2(1, 1) = — e*'^/^ and ap {j, k) = I otherwise is a set of orthogonal one dimensional 
projection operators |l9|. If we make this definition for d not prime, we find that we generate 
rank 1 projection operators which are not orthogonal. The reason that the factor 0:2 appears 
in the p = 2 case is that for p an odd prime S'J^ = 5*0,0, however, (02(1, 1)5*1,1)^ = 5*0,0 since 
Sfi = — 5o,o- 
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IV. MUTUALLY UNBIASED BASES I 



We review the theory of mutually unbiased bases (MUB) for a particle whose state vectors 
lie in a p-dimensional complex Hilbert space Hp, where p is a prime. It can be shown that 
there exist p + 1 orthonomal bases (ONB) in this space which are MUB |ll|,Q,Q; that 
is, if ip and are state vectors that belong to a pair of ONB that are mutually unbiased, 
then I (01?/') I = The simplest example of mutually unbiased bases occurs for p = 2, for 

which the bases are composed of the eigenvectors of the three Pauli matrices {crx,<Jy,<Jz}- 

There is a nice way to characterize the MUB using commuting classes of the generalized 
spin matrices This leads to a natural way to introduce discrete phase space, and, in 
turn, to a definition of a Wigner function. We denote the two dimensional vector space with 
components in Zp by V2{p), and use the letters u and v to denote vectors in this space. This 
vector space contains p"^ distinct points, and it is convenient to index the spin matrices 
using V2{p), 

V = {VQ,Vi) ^ Sy = Syg^y^. (16) 

With this notation Eq. (fT^ becomes 

SuSySi = Ti^^'^Sy. (17) 

It follows from this that two spin matrices commute if and only if the symplectic inner 
product of their index vectors vanish. Therefore, the problem of finding commuting sets of 
operators is transformed into finding solutions to the equation u o v = Q for vectors in the 
two dimensional vector space V2{p)- The solutions are easy to find; the p + 1 index vectors 
Ma, a G /p = {0, 1, ■ ■ ■ ,p} partition the spin matrices into p + I sets defined by 

Ca = {bua = 6(1, a), beZp}^ma = {S^, beZp} a<p 

Cp = {hup = 6(0, 1), beZp}^mp = {Si, 6 G Zp}. (18) 



(Note: in j20| Cp was denoted by Coo). 

Equation p8|l relates each vector in V2{p) to commuting sets of unitary matrices such 
that dJta n DJlb = {>S'o,o} for a 7^ 6 . This follows from the fact that in V2{p) two non-zero 
vectors with vanishing symplectic product must be coUinear. The state vectors in each basis 
are the eigenvectors of the associated set of unitary matrices in Eq. (|18|). The projection 
operators for these vectors are defined in Eq. (fT^ and can be found in ll|. 



► 
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FIG. 1: The vectors Ua in V2(3). 

V2(p) will be used as the phase space for a single system with Hilbert space Hp, and vectors 
in V2{p) will be used as indices for the characteristic function and for the Wigner function. 
The "horizontal" and "vertical" axes of V2{p) are associated with the spin matrices Suq 
and Supi respectively. In general, a vector (or point) (j, k) in V2{p) corresponds to 5*^^^. The 
projectors generated by Suq are associated with the basis is j G Zp], and the projectors 
generated by Su^ are associated with the basis {\k) = ( ij A ^p) J^ ^j-n v'^-^li), k G Zp}. The 
latter states are often referred to as the phase states, j2c 



27| . The Hermitian operators 
J = X]j=oili)OI with eigenstates {|j)} and $ = J2^=ok\k){k\ with eigenstates {|A;)} are 
said to be conjugate observables, since these states are Fourier transforms of one another. 
This is in analogy with the operators Q and P of section IH] although the commutation 
relation of J and $ is not proportional to the identity operator, and is, therefore, state 
dependent. 

The fact that the sets Ca correspond to a set of MUB can be seen by computing the 
projection operators for the sets, and showing that [2^ 

$^P«„(r) = ^0,0 (19) 

tr[P^Sr)PuSs)] = 5{r,s) (20) 
tr[P^Sr)PuSs)] = ^ for a ^6. (21) 
In particular, the proof of Eq. (pT| depends on the orthogonality of the spin matrices and 
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the fact that tr{Sj^k) = for all the spin matrices except the identity. This set of MUB is 
complete in the sense that there are p + 1 ONB in the set, the maximum number possible 

Q. 

V. THE DISCRETE WIGNER FUNCTION FOR A SINGLE PARTICLE 
A. The Wigner function 

The discrete Wigner function of interest here was introduced by Wootters in Fol- 
lowing Wootters we wish to define the discrete analog of the Wigner function such that 
properties 1 — 5 of section [H] are preserved. Our approach differs by emphasizing the role of 
the spin matrices. 

Let p be a prime number, and p G Mp be a density matrix describing the state of a system 
on the Hilbert space Hp. Define the characteristic function over V2{p) 

Xpimua) = ximua) = tr [p {ap{ua)SuJ'^] , (22) 

where ap{u) is defined above Eq. p5|) . The properties of x that we shall need are 

X(0) = 1 (23) 
X{mua)* = x(-mMa) (24) 

This last result follows from the fact that (S*^)^ = S~"^, since Su is unitary. 

Let V = {vq,Vi) and u = {uo,Ui) be vectors in V2{p). Then using Eq. (fTTj) the discrete 
Wigner function is defined as the discrete symplectic Fourier transform of the characteristic 
function: 

Wp{v) = W{v) = 4 E ^""^(") 

= -2 U(0) + E E ^™"x(mn,) ] . (25) 

P \ a=0 m=l / 

The sum over m excludes the m = term, which gives rise to the first term in brackets. 
The equality of these two expressions follows from the fact that the vectors {mua, a G Ip, 
m G Z*} U {(0, 0)} = V2{p), where Z* = Zp — {0}, that is, these vectors partition the space 
into distinct lines through the origin. This fact illustrates the role of the geometry of V2{p), 
see appendix IX El 



If we substitute Eq. into and use (fT3jl. the Wigner function can also be written 



as 



[voUa)] 



(26) 



where {pr(r\a,p) = Tr [p-Pu„(r)] , r G Zp} is the probabihty distribution that can be esti- 
mated from one of the p+1 experiments determined by the set of MUB j2^. The sum over a 
gives a complete set of measurements for determining the Wigner function or, equivalently, 
as we shall see, the density matrix. This form of W shows that it is real and that it may be 
negative. 

Equation can be rewritten as 



W{v) 
A{v) 



tr [pA{v)] 



V O U„ 



0,0 



(27) 
to define 



The set of Hermitian operators {A{y), v G V2{jP')} was used by Wootters in 
the Wigner function and is an orthogonal basis of Mp. To verify this, one uses the MUB 
properties from Eq. (pT| and computes as in 



Tr [A {u) A {v)] = \ 



V 



2 (p + 1) + J] J] Tr [P„„ {u o ua) Pu, {v o u,)] 



Regardless of u and v each term in the double sum equals 1/p when a ^ h. \i a = h and 
u = V, each of the resulting p+1 terms equals 1. If a = 6 and u ^ v, then the trace equals 
zero except for the one case when u — v = mUa so that uoUa = v oUa and the trace equals 
one. Collecting terms gives 



Tr \A{u)A{v))] = -5{u,v). 

P 

Note, by the way, that one can use the orthogonality to express the identity as 



(28) 



(29) 



In the preceding discussion we have written the Wigner function and the characteristic 
function. In fact, for a given density matrix and a complete set of MUB, a class of Wigner 
and characteristic functions can be defined. For example we can multiply the characteristic 
funtion in Eq. (j22p by an appropriate phase factor and get a new characteristic function 

Xp{mUa) Xp,ra{mUa) = T]'^'''' Xpi^Ua) , 
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where G Z^. Under this transformation 



1 / " 



where r = (rg, . . . ,rp). This approach provides an operational way of defining the class of 
Wigner functions described in 3] and in the recent work of 0]. 

Before showing that the definition Eq. has the desired properties, we present three 
examples. 

B. Examples 

1. Qubits (p=2) 

Using Eq. (jHl), the spin matrices may be shown to be equivalent to the Pauli matrices: 

where (Tq is the 2x2 identity. The classes of MUB are generated by 

Co = {6(1,0)} ^{ao,aJ 
Ci = {6(l,l)}-.{oro,«a,} 
C2 = {6(0,1)} ^{ao,f7,.} 

where 6 G ^'2. The most general density matrix may be written as 



P 2 

where {mx,my,mz) is a vector with real components and length less than or equal to 1. In 
this case 

Xiuo) = m^, xiui) = my, ^(^2) = m^- 

We have included the factor a2{u) so that x is real. For p = 2 we have rj = —1, and for 
V = {vo,Vi) G V2{2) 

W{v) = ^ (1 + ^zV""' + myT]^"^-'"'') + m^.r^-'^o) . 
11 



It is now easy to see that summing over a horizontal hne gives 



1 



PuM 



where -P?io(0) is the projection operator for the state polarized along the positive z-axis, and 
Puo(l) is the projection for the state polarized along the negative 2;-axis. A similar result 
holds for the sum over a vertical line, that is, a sum over vi and the x -axis. For s G Z2, 



which corresponds to summing along the line 1), h G Z2}. Finally, for this case, the 
Hermitian matrices defined in Eq. ()27|1 are 



It is well-known that for a single particle W{v) can serve as a hidden variable probablility 
distribution if it is nonnegative. This is because, as we shall see below, the measurement of 
an arbitrary observable O is given by 



where Wo{v) is the Wigner function defined with p replaced by O in Eq. (j22I)- Therefore, if 
the Wigner function is non-negative we can construct a complete hidden variable theory of 
a single qubit consistent with quantum mechanics. However, there appear to be cases where 
this does not work, for our present example if m = (1, 1, l)/-\/3, then 1^(0,0) < 0], see 
however j3] where it is shown that a hidden variable theory can always be constructed for 
a single spin. We also note that since this m corresponds to a pure state there are bases 
in which W{v) > 0. The positivity of the Wigner function is therefore sufficient but not 
necessary for the existence of a hidden variable theory. For a discussion of the positivity of 
the Wigner function see Q]- 



V 





tr{pO)=p W{v)Wo{v), 



feV2(p) 
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2. A pure state in Hp (p > 2) 



Let p — Pu^{r), then for a & Ip 



i|:,-"tr(sr.(sy) 

^ k=0 



1 



S{a,b)r) 



Therefore, 



W{v) 




m=l 



,li;0-MaOl'l) 



1 



S{r + UaOv,0) , 



P 



so that W{v) vanishes except at points along a hne in V2{p). In particular for the case a — p, 
W{v) vanishes everywhere except along the vertical line vq = constant, and for a = 0, W{v) 
is constant along the horizontal line Vi = constant and vanishes everywhere else. 

Given an arbitrary pure state, we can always find a MUB that contains this state as one 
of the basis vectors. This shows that there is always a MUB for which a pure state has a 
non-negative Wigner function. On the other hand if the pure state is not chosen as one of 
the MUB vectors the result is more comphcated as will be seen in example 4 below. 

3. Completely random state 

The density matrix for the completely random state is p = which gives W{v) = 

l/p^i that is, W{y) is constant. 

4- The operator O = \ j){k\ 

As stated above, we can define a Wigner function for operators other than density matri- 
ces. We give an example here which we shall use later. For the case that p is an odd prime, 
let \j) and \k) be vectors in the standard basis and let 



0^\j){k\. 
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Then for a < p 

Xj,k {mua) = Tr [\j) {k\ {SuJ""] , 
and, for reasons that are explained in Section I VI II we introduce a phase factor when a=p 



r / -1 \ ' 
Xj,k (mup) = Tr {k\ [t]'^ Su^j 



where —2 ^ is taken as {p — l)/2 since in the exponent we can compute modp. Using 
Eqs. dH) and (II2D, we find 

Xj,k (mup) = rf '('^^^^^ (j, m + k). 

Working through the details gives 

WuKki (v) = (0, v^ + (j + k)) , (30) 

where v = (fo,fi). Note that if = j, is a density and is a special case of 

example 2 above. For j ^ k, we get 



Woiv)* = Wo^{v). 
Now suppose that = J2^ZoCj\j), then 

1 



P r=0 

As stated above this is a more complicated form than we found for the case = Pb{r). 

For the case p = 3 we illustrate this in Fig. 2 for the case \ip) = + |2)). 

C. Properties of the Discrete Wigner Function 

We now examine whether the definition ()25p or, equivalently, ()26p satisfies the criteria 
that we set out in part 1. 

1. The mapping is p ^ Wp is linear on and convex linear on the density matrices. 

2. W{v) is normalized since 

p-i ^ / p p-i \ 

J2 ^(^) = -[p' + J2Yl ximUa)p'S{m, 0) = 1. 

vovi=0 P \ 1=0 m=l / 
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FIG. 2: The Wigner function for p = 3 for the pure state |^) = + |2)). 

3. The reality of W follows immediately from Eq. (f^Bjl . 

The first three results also follow directly from Eqs. (f?fj) and ((201) • 

4. For w G V2{p), if p' = Sl^pSw then using Eq. (fTTj) 

Wpi{v) = Wp{v + w). 

Note that if p commutes with S^, the Wigner function is invariant under translations along 
w. Furthermore, the characteristic function vanishes for u such that w o u mod p. 

5. The marginal distributions are easily computed. We consider the more general case 
of summing along the points on any of the lines in phase space, where a line in phase space 
V2{p) is defined as the set of points that satisfy the equation 

L{b, s) = {{x, y) : -xh + y - s = {x,y) o Ub - s = 0, b,s e Zp} 

L{p, s) = {{s,y) : -X + s = {x, y) oUp + s = 0, b,s e Zp}. (31) 

L{b, s) is the line with "slope" b which intersects the vertical axis at s, and L{p, s) is a 
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"vertical" line that intersects the horizontal axis at x = s (see Appendix IX E|) . Let 

Ms) = Y,W{v)6{vou,-s,0), 



then using Eqs. (|26p and (|2ip. we can show that 



tr 



r=0 



+ ptr [pPu,{s)] 



tr[pPu,{s)] 



We have used the fact that for a ^ b the sum over v becomes a sum over Zp, and this sum 
is the identity operator, while for a = b, we have voub = s. Therefore, we see that summing 
the Wigner function over any line in phase space gives the probability that the system is in 
the corresponding MUB state. 

6. Since Wp and xp s^i'fi Fourier transforms of one another, Plancheral's formula gives 

p p-i 



P' 



(32) 



a=0 m=l 



We also have, setting mua = (j, k) = v, 



\Xpiv)\^ = \trpS^\'^ = {p,S^){Sy,p). 



using Eq. ((7j). More generally, 

Xpi(^)Xp2(^) = (trpiSy)* trp2S^ = {pi, S,a){Sy, P2). 
Summing over the complete set of S^, from Eq. (^U, we can write Plancheral's formula. 



tr [piP2]=pJ2Wp,{v)Wp,{v). 

V 

See also js| where the derivation is based on Eq. (j2HI)- 

The support of a function /(t>) on phase space is defined by 



(33) 



suppif) = {vE V2{p) : f{v)^ 0} 



(34) 



and \supp{f)\ is defined as the number of points in supp{f). From Eq. we have 

pW;ivo)<pJ2^pi^)=^rp^ < 1- 



(35) 
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which imphes that for any point Vq, 
Now let 

{1 V E suppWp 
otherwise. 

Then applying the Schwarz inequality to the normalization equation and using Eq. ()35p we 
get 

1 = J2Wp{v)nw,{v) < /Xl^p^^)'"^p(^)Zl^^p(") - J-\suppWp\ 

V ^ V u \ P 

or 

\suppWp\ > p. 

This is analogous to the continuous case where the uncertainty principle implies that W{q,p) 
can not be concentrated into too small a region. We have seen that if p is a pure state selected 
from the MUB that \suppWp\ = p, and Wp{v) = 1/p on its support, so the lower bound is 
attained. 

If pi and p2 correspond to orthogonal states, then Eq. gives 



J2Wp,{v)Wp,{v) = 0, 



which along with the normalization condition implies that, either suppWp^ and suppWp^ are 
disjoint or at least one of the Wigner functions must take on negative values. For example, 
we saw in IVB 2l that the orthogonal states in one of the bases of a set of MUB have support 
on non-intersecting lines of V2(p). 

There is an inequality, referred to as an uncertainty principle, that also follows from the 
discrete Fourier transform: 

IswppVFpl \suppxp\ > P^, 



22| | . Equality holds for the random state discussed in example IVB 31 above. 



D. Inversion formula 



In the case of continuous phase space, the density matrix for a particle confined to one- 
dimension can be obtained from Eq. ^ by using the inverse Fourier integral. We can 
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proceed in a similiar manner for the discrete case. First using the discrete Fourier inversion 
formula, 

X{mua)= W{v)r~'^ (36) 

Then from Eq. (j22j) and the completeness of the spin matrices 

P = I ( ^0-0 + E E ^ ("^"")* ^uA ■ (37) 
Substituting pUI) into P7j). and using Eqs. (fT^ and (|77j). we also get 

p= 5^ W{v)A{v). (38) 

Therefore, we have an expression for the density matrix as an expansion in the spin matrices 
with coefficients given by the characteristic function and an equivalent expansion in terms 
of a basis of Hermitian operators with the Wigner function as coefficients. 



VI. MUTUALLY UNBIASED BASES II 



To define the generalized spin matrices in the case when c? = where p is prime, 
we require the notion of a finite or Galois field GFlp"-), see Appendices IX Al and IX Bl for 
more details. There is a systematic way of representing the elements in GF (p") that uses 
the structure of polynomials irreducible over GF {p) = Zp. An irreducible polynomial is a 
polynomial f{x) of degree n with coefficients in GF{p) that can not be factored into non- 
constant polynomials of lower degree. Then the elements of GFijf-) may be represented by 
polynomials of degree less than n with coefficients in GF{p). The simplest example is that 
of two qubits, p = 2, n = 2. In this case the irreducible polynomial is unique and is given 
by + a; + 1. Define A to be a symbolic solution of + x + 1 = mod2. Then every 
element of GF{2'^) can be written as a = Oq + Aai where and ai are in GF{2) . This 
is analogous when working with real numbers to letting i denote a symbolic solution of the 
equation x^ + 1 = and introducing complex numbers as x + iy. 

For the case of n = 2 and p an odd prime, let D be an element in GF {p^) such that there 
is no solution in Zp = GF{p) to the equation x^ — D = modp. In technical terms, D is a 
quadratic non-residue of p. There are an equal number of quadratic residues and quadratic 
non- residues in GF[p). Then elements in GF (p^) can be represented as j + kX, where j and 
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k are in GF (p) and A is taken to be a symbolic solution of — D = modp . Addition 
and multiplication of elements of GF{p^) are defined by 



(ji + fciA) + (j2 + fcaA) = (ji+j2) + (fci + fc2)A 

(ji + /i;iA)(j2 + /i;2A) = {3132 + Dkik2) + Uik2 + ki]2)\ 

where the additions in the parentheses are modulo p. We refer to Appendix IX Bl for more 
details. 

We can construct a complete set of mutually unbiased bases when d = by following 
the same procedure that was used in the d = p case |i2B]. The key idea for constructing a 
MUB is based on the fact that we can define a two-dimensional vector space V2(p") over 
GFIp"-), and + 1 generating vectors Ua where a is in the index set Ipn = GFijf') U {p"}. 
Specifically, define 

/ (l,a), aeGF{p-) 
Ua= { (39) 

y (0,1), « = p" 

Each of these vectors can be used to define a class containing p" vectors, 

C^ = {l3u^,(3eGF{p^)} (40) 

where a G Ip-n. Each pair of vectors in a class has vanishing symplectic product, Eq. (|lip 
where the operations are with respect to GF{p'^). We want to find a spin matrix representa- 
tion of these classes, that is, we wish to find a mapping from this space to the set of tensor 
products 

n-1 

5u = (g)5,. (41) 



r=0 

n-1 



where u = m*-'^'-' G V2n{p) and each u^"^^ G V2 ip) = V2{p). To do this we define an 

r=0 

isomorphism 

n-l 

M:V2{pn^V2n{p) = ^Vi'\p) 

j=0 

that preserves the symplectic product in the following sense. For each vector v G V2{p^), if 
M{v) = V = @^ZqU^-'\ where u^^^ eV2''\p) , define 

p-i 

^r,o^r2 = Y,i^? °^2^)- (42) 

j=0 
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Then f i o t>2 = implies vi o V2 = as is shown in Eq. ()87|) . We present an outhne of the 
derivation of M in Appendix IX Bl and refer to [20] for another discussion. It is worth noting 
that the construction of V2n{p) is analogous to what is done in the continuous case. There 
we take the direct sum of the two-dimensional vector spaces corresponding to independent 
conjugate position and momentum pairs. 

To perform the analog of what was done in Eq. ()18|) . it is useful to introduce the generators 
of the index set for the set of MUB, again the details are given in Appendix IX Bl For 
define X'Ua G V2(p"), r = 0, ■ ■ ■ , n — 1. Then define a set of generators on V2„(p) 

= {gr(«) = M {yu^) , r = 0, ■ ■ ■ , n - 1}, (43) 

and define the corresponding spin matrix using Eq. (PT| as 

n— 1 

^i(«)-(8)^:<^.(.) (44) 

i=o 

where each S^u) acts on a Hilbert space Hp. The generalization of Eq. (fTHj) is 

{n— 1 n— 1 n— 1 ^ 

n^t(.)=(8)n^!o)(.)' ^reGF{p)\ (45) 
r=0 j=0 r=0 ) 

for the generation of + 1 disjoint sets of of commuting operators 93Tq, where OJIq, flOJl/j = 
{5*0.0} for all a ^ (3. We have written the mapping in Eq. (jl^ from the set of basis vectors 
Ga rather than the space C^. 

It is also possible to write down the set {Pa {r) ,r ^ Vn (p)} of rank one orthogonal pro- 
jections defined by each of the + 1 commuting classes 3JIq,. This gives the set of MUB 
as projections defined explicitly in terms of sums of the spin matrices in each class. The 
procedure to do this is discussed in j^], and is illustrated there for the case for n = 2 . The 
corresponding projection operators for the case p > 2 are 



r=0 

Pm^M) = (46) 

br=0 J 

where s = (sq, . . . , s„_i). For p = 2 it is necessary to include the factors a2{j,k) in the 
definition of the projection operators as shown in Eq. (|15p. (r) has trace one, and it is 
straightforward to check that if r 7^ s 

P4r)P„(s) = (5(r,s)P4r). 
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It is easy to show that each Pa (r) is a product of commuting projections. One can also 
show that Pa (r) = {Pa and it follows that Pa (r) is a rank one orthogonal projection 
and that 

r 

Finally, it can be shown that for a j3 that 



tr [Pa{r)Pp{s)] = 




The explicit calculation of the projections and of the set of vectors in {p) corresponding 
to Ca depends quite specifically on p and n and on the representation of elements in the 
different finite fields. When n = 2 and p is an odd prime, however, one can give a unified 
summary of the results of the theory. Without going through the detailed construction 
outlined in Appendix IX B[ it is easy to check that the vectors in each of the classes below 
have symplectic product zero. 
Example d=p^ 

Let d = p^ with p an odd prime and D such that — D = Q modp has no solution in 
GF{p). In Appendix IX Bl Eq. (j^^ it is shown that the p^ + 1 commuting classes of indices 
are generated by 

Ga,,a, = {(1, 2ao, 0, 2Dai) , (0, 2Da^, 1, 2Dao))} 
where ao and ai are in GF{p), and 

Gp2 = {(0,1, 0,0), (0,0, 0,1))}. 

One can check directly that the vectors in each Ga^ai have vanishing symplectic product. 
Then the spin matrices that generate the commuting classes may be written as 

\ {Sl,2ao ® So,2Dai) " (5'o,2Dai ® Si^2Dao) ^ , ^0, h G GF{p) > , 

Gp2 ^ mp2 = { {S',\ ® 5-0,0) (5-0,0 ® 5'o':i) ,bo,he GF{p) ] . (47) 
The corresponding projections are given by 

Pao,a. ((ro, ri)) = -Y, {iv''S,,2ao ® So,2Da,)''') ' {{v''So,2Da, ® S,,,naof') 

^ bo=0 ^ bi=0 

Pp- ((ro, ri)) = -Y, {iv^'So,i ® 5o,o)'°) - ((^''^o-o ® " 

^ 60=0 ^ bl=0 
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We note that each of these one-dimensional projection operators is the product of two 
commuting rank p-dimenional projections. The two p -dimensional spaces that they project 
onto intersect in a one-dimensional space. 

VII. WIGNER FUNCTION FOR d = p'^ 

In earlier work P, IsJ , the phase space on which the Wigner functions were defined when 
d = was chosen to be V2(p"). The advantage of this choice is that one can use the 
underlying geometry to great advantage. The disadvantage is that one has to label coordi- 
nates using elements from the Galois field GF{p^) which does not lend itself to a discussion 
of separability. However, as we saw in Section IVIl and as is elaborated in Appendix IX Bt 
there is a natural isomorphism M between V2(p") and V2n{p) which encodes the geometry 
of V2(p") in V2n{p)- We take advantage of this structure to define our Wigner function on 
V2n{p) This is in close analogy to the continuous case and simplifies computations involving 
the generalized spin matrices. 

In particular, this approach enables questions involving separability to be treated effi- 
ciently. In this section we illustrate the ideas in detail for n = 2, leaving the generalizations 
to the next section and the Appendix. 

A. Separability of the Wigner Function for p an odd prime 

We consider a bipartite system composed of subsystems of dimension p, a prime. As we 
saw in Sectional there is a certain latitude in the definition of the Wigner function that is 
available because of the freedom to include phase factors in the characteristic function. Our 
goal in this section is to show how that freedom enables us to define Wigner functions for 
one and two subsystems so that separability is respected. Specifically, for a product state 
we want 

p = r(^^^Wp{u) = Wr (?/(°^) Wp , (48) 

where u = u^^^ © u^-^\ Then, since Wp (u) is convex linear on the space of densities, we will 
have the general statement that 

P = Y^PkTk ®Pk^Wp (u) = Y^P^^r^ (^^°^) («^'^) • (49) 

k k 
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A natural definition of tlie cfiaracteristic function x = Xp is to use Eq. 
and define 



witli n 



tr 



(a) 



(50) 



where w = 6ogo(«) + bigi{a). We can rewrite tlie product of the S matrices on Hp2 as a 
direct product of S matrices on Hp^^ ® Hp'^\ Rewriting w = ■u*^'^^ © u^^^ where ■u*^-'^ G V2''\p), 
we find 

X(w) = r]^tr [pS^m ® S^m] . 

The problem with this definition is that in general $ 7^ 0('U*^°)) + </)('u'^^^), so that the corre- 
sponding Wigner function would not factor when p is a product state. Now as pointed out 
before, there is some freedom in the choice of phase in defining the characteristic function 
and the Wigner function. For this reason it is convenient to introduce a phase factor into 
the definition of the characteristic function to avoid this problem. We shall therefore define 
the characteristic function as 

X(w) = r/-®x(w), (51) 

using Eq. and the B defined in Eq. that is linear in 60 and 61. The linearity in the 
b's is important, as we shall see, because we want to write the analog of Eq. with the 
appropriate projection operators given in Eq. 

The underlying reason for having to introduce the phases arises from the fact that we 
are using the geometries of ^(p^) and V4(p). That fact forces us to go into some detail to 
define appropriate phase factors and to confirm that they work. 

For example, consider the case of p odd discussed at the end of the last section. For 
a ^ define 



w 



Tr 



p(r^-^-5,„)^°(r^-5,0^^ 



(52) 



and for a = p^ define 



Then define 



Xp (w) = Tr 



(53) 



(54) 
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Note again that — is computed modulo p and equals {p — 1) /2. The vector symplectic 
product in the exponent of 1] is defined in Eq. (021) • From Eq. (^3)) we can write out the 
right hand side of Eq. (j48p with one modification. For mup = m (0, 1) take 

(55) 



Xrimup) = tr 1^7] Su^ 

as we did in Example 4 of Section^ Then Wr (u) is defined as the usual symplectic tranform 
and can be written as 



Wr (U) = - 



p—1 p^l p^l 

.m=0 c=0 m=l 



(56) 



Finally, we get the right hand side of Eq. pHjl as the trace oi {j ® jj) times the expression 

p-i p-i 

mo=0 mi=0 

^uo(mo«co©miMci)^-2-i(5(co,p)mo+5(ci,p)mi) )'"°®(5'« )'"^. 

(co,ci)^(p,p) (mo,mi)^(0,0) 

The left hand side of Eq. (pSj) can be written as the trace of (r (g) /i) times the expression 



p— 1 p— 1 
60=061=0 



^ ^ ^uo(6ogoW+6igi(a))^-Dai(6o+^>i) (5g„(„))'" (5g,(,))''^ . 
«7^p2 (bo,fei)7^(0,0) 

Note that in this equation we have the ordinary matrix product in the second term. 

Our goal is to confirm that Eq. (j48j) holds with the above definitions of the characteristic 
functions. Using Eq. (jl2p we can pair the indices of the spin matrices in Eq. (j48|) to obtain 
the index equation relating terms in rhs to Ihs, 



w = moUco © miMci = &ogo (a) + &igi (a) 



(57) 



which includes the w = (0, 0, 0, 0) term that is incorporated in the first summations. It 
follows that the phase factor 7^"°(^ogo(a!)+f'igi(a)) jg common to the corresponding terms of 
rhs and Ihs, and we can cancel it. It is also obvious that the a = p^ terms equal the 
corresponding terms associated with cq = Ci = p and that the remaining phase factors in 
this case are also equal if we set = 6^. 
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To match terms in the second sets of summations, we multiply out the powers of the spin 
matrices in Ihs to obtain 



SmoUcQ ® SmiUc-^ — Sbo(l,2ao)+bi{0,2Dai) ® 'S'feQ(o,2Dai)+fei (l,2Dao) ) 

where the equality follows from the index equation. This process introduces phase factors 
using Eqs.()10p and (fT^ . and it remains to prove that the resulting exponents of t] are equal. 
Specifically, one has to verify that subject to Eq. 



(1 -5(co,p))co( 2 1 +(l-5(ci,p))ci( ) -2 ^(5(co,p)mo + 5(ci,p)mi) (58) 



equals 



-Dai {bo + bi) + 2ao ( ^ ) + 2^«o (2 ) + ^b^biDai. (59) 



We verify the equality for a 7^ by considering different cases. Let a = aQ + ai\. If 60 and 
61 are both non-zero, tjiq = bo, mi = bi and 

boCo = bo2aQ + bi2Dai, 
bici = {bo2Dai + bi2Dao) ■ 

If bo = 0, 

w = {0,bi2Dai) © {bi,bi2Dao) 
and Co = p, mo = bi2Dai, ci = 2Dao and mi = 61. Similarly, for 61 = 

w= {bo,bo2ao)®iO,bo2Dai), 

and mo = 2Dao, m^ = bo, m^ = p and mi = bo2Dai. Substituting these expressions in 
Eq. gives (j^ . We have gone through this in some detail because the method illustrated 
generalizes to the case of complete separability of n subsystems. It should be noted that the 
argument leading to Eq. ()48p did not require that r or /i be a density matrix. 

Our ability to add a phase factor to the definition of the characteristic function is related 
to an arbitrariness in the assigning of state vectors in a basis on the Hilbert space to lines 
in phase space as noted in ^j. This is illustrated in I VII C 2) below. 

A different definition of the Wigner function in terms of the characteristic function can 



be found in 



2jj . Vourdas replaces the M transformation by introducing the trace operation 



into the Fourier transformation. 
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B. Properties of the Wigner Function 



Because we have used the same format in defining the Wigner function for two subsystems, 
Eq. (j^H), as was used in defining it for a single subsystem, Eq. we expect the properties 
in Section ini to hold. With the definition of Xp(w) in Eqs. and (fHH|) conditions (PHj) 
and are satisfied. The discrete Wigner function Wp for a density p on Hp2 is defined 
using the symplectic Fourier transform Eq. (jSH); consequently, Wp is convex linear on the 
space of densities and linear on the space of p"^ x p"^ matrices. Again, the defining Eq. ()54|) 
is invertible, so that one can obtain the Xp (w) and thus the spin coefficients of p from the 
Wigner function. With this definition Plancheral's formula becomes 

veViip) q,r=0 M 

We also have, as in Eq. (jH^ . that 



trM=p'J2Wp,MWp,M 



and, consequently, |Vr(v)| < 1/p and [^^^^^^(v) | > Using the notation of Eq. (UBI), 

we can write 

Wp{u)=tr[pA{u)] (60) 

where 

p^A (u) = -^0,0 ® ^0,0 + Pp^ (-2-' + u o go , -2-1 + u o (p^)) (61) 
+ E Pc. {-Dai + u o go (a) , -Dai + u o g^ (a)) , 

corresponding to Eq. ()27|) . From Eq. ()61|) it follows that VFp is real for densities p. In 
particular, {y4(u)} again defines a complete orthogonal set of Hermitian matrices. The 
argument is analogous to that leading to Eq. (27) and leads to 

Tr [A (u) A (v)] = p-^6 (u, v) 

Thus we can interpret the Wigner function Wp as the set of coefficients of p in the orthogonal 
expansion relative to {yl(u)} analogous to Eq. (jHHj) . 

The analogues of the other properties of Section IV CI follow in the same way as before. 
W is normalized since we can use Eq. fl60|) to prove 



J2Wp{u)=Tr p5^A(u) =Tr[p] 
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since Eq. holds in this case. If p' = SlpSz E Hp® Hp, then Xp' (w) = 'rf°^xp (w), and 
Wp' (u) = Wp{\i + z) as before. 

Summing Wp (u) over a "hne" in V2 {p^) corresponds to summing over a translation of a 
two dimensional subspace in V4 {p) and again leads to a marginal probability Tr [pPa (si, S2)]. 
To see this let Ca denote the two dimensional subspace associated with a. It is easy to show 
that 



This can be written as the trace of p against the projection Pq(si,S2) for appropriate 
indices si and S2 which depend on r and the phase factors used in the definition of the 
characteristic functions. Thus using definition of Eq. the Wigner function satisfies the 
conditions proved in section IVT^ and the requirement that Wp factor for separable p as in 



As pointed out to us by Wootters, Eq. may be used to give a positive answer to a 
question posed in j^. That is, with the phase factors given above, we have 

where u = u^'^^ © m*-^-*. The proof is easy, rewrite Eq. ()48|) as 

tr [to ® riA(u)] = tr [tq ® ^^(m^") © A{u^'^^)\ . 

This equality holds even the r's are not densities. Since Hermitian matrices of the form 
T ® p form a basis of Mp2, this inequality holds for all ^(u) for all u G Vi{p). 

C. Examples 

1. Maximally entangled state 

For prime p let |^) = ^ |j) Ij), so that p = |^) (^| = 1 Y.j,k lj) (^1 ® lj) (^1- % the 
separability property and linearity we know that if u = u^^^ © u^^"* = (xq, yo, xi, yi), then 




Eq. (gni). 




where W^|j)(A:| (u) is defined in Eq. (jHIUl . It follows that 



Wp H = ^Y1 {yo + + k),0)6 [yo + 2'' (j + A;) , O) 
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and simplifying we get 

Wp (u) = ^5 (1 + xo + xi, 0) 5 (yo, yi) 

Thus the Wigner function for this maximally entangled state equals l/p^ for the p"^ four- 
vectors with u^^^ = {xo,yo) and u^-^^ = (—1 — xo,?/o) and equals zero elsewhere. 

Although the Wigner function for this state is positive, it is a non-classical state. In 
particular, entangled states violate Bell inequalities. Since the Wigner function discussed 
in this example is not separable, it need not respect mathematical inequalities based on 
separability. 



2. MUB 



Let p = Pa {sq, Si) . In this case it is simplest to use Eq. flUTj) so that 
Wp{u) = ^tr[Pa{so,Si)A{u)] 



P 



-l + ^Tr[Pa {sq, Si) Pfs (ro,/? (u) , ri^fs (u))] + Tr [P^ (sq, Si) Pa (ro,„ (u) , ri,„ (u))] 



= ~2 [""^ + pVp^ + ^ (■50, ro,a (u)) 6 {si,ri^a (u))] 
= ^6 {so,ro^a{u)) 6 {si,ri^aiu)) 

Thus, Wp{u) equals on those four- vectors which match the given phases and 
equals zero elsewhere. For a ^ p^, r\^^a (u) = —Dai + u o (a) , and it can be shown easily 
that the set of four-vectors satisfying those conditions is 

{u = bogo (a) + 6igi (a) + (0, Sq + Dai, 0, Si + Dai) ■ &o, h G GF{p)} . 

That is, Wp (u) is constant on a shift of the two-dimensional subspace indexed by a. An 
analogous result holds if a = p^, and, as expected, this parallels the situation when n = 1. 



D. Separability of the Wigner function for p=2 

When p = 2 Eq. ()54|) can be used to define the Wigner function with the definition of the 
characteristic function given in Eq. (jU^ below. Properties other than separability follow as 
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before, but the analysis leading to separability for j9 odd does not work in this case. The 
discussion above made use of the existence of a quadratic non-residue D\, however, for p = 2 
no such quantity exists. In addition we must include the factors of ^2 = Q!2(l, 1) = —i 
defined at the end of 11111 

Explicit forms of generating vectors are 

Gao,ai = {(1, ai, 0, ao + ai), (0, oq + ai, 1, ao)} 

for a = aQ + aiX E GF{2'^), and 

G4 = {(1,0, 0,0), (0,0, 1,0)}. 

For the case a 7^ 2^, the analog of Eq. (j^H) is 

Xp(w) = Tr[p(arr^'^"5gJ^««r/^'^5gJ^^] (62) 

where tq and ri depend on ao and ai. It is convenient to write the index equation Eq. (jKTjl 
in the form 

w = (60, go) © (bi, qi) = feogo(a) + &igi(a), 
then it is not difficult to show that for {bo, bi) 7^ (0, 0) 

ao = boqo + (60 + &i)gi 
fli = (bo + bi)qo + biQi. 

This allows us to replace the sums in the Wigner function over and ai by sums over go 
and qi. Now we can write 

As stated above, we require that the phase factors are linear in the b's. In order to enforce 
this it is easy to show that if ro = and ri = = boqo + &i(cio + ^i) the exponent of f] is 
simply 61 gi. This calculation makes use of the binary arithmetic, in particular 6^ = b. 

Finally, we find that for p = t ^ fi the phase factor t]^^'^^ requires that we use different 
one particle Wigner functions for the two particles. Equivalently, 

where /z* is the transpose of the qubit density matrix /i. If we had taken tq = ai and ri = 
the transpose would have appeared on r, rather than on p. 
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E. Separability and Partial Transposition 

A nece.a.y condition fo. separability of a density n.at.ix of a bipartite system p . H^^H, 
is the Peres condition [18||. That is, the density matrix must transform into a density matrix 
under partial transpose 

PT : {jo,3i\p\ko,h)^{jo,h\p\koJi) (64) 

The transpose of a spin matrix is given by {Sj^kY = V~''''Sj,p-k', consequently, under the 
PT transformation 

Therefore, 

PT : Ty(u) =Ty ^ TV - (n« + 1)) 

Unfortunately, this is not very useful since proving that W corresponds to a density matrix 
is not simple, see IX Gl 

VIII. WIGNER FUNCTION: d = p^. 

The generalization to degrees of freedom, where p is prime, is based on the Galois field 
(see GF{p"') and Appendix IX Bj) . Starting from Eqs. and pn|l . the set of vectors 
in Ca defined on the phase space V2(p") generates a MUB. As before u denotes a vector 
in (p) = © (p) that we also write as u = ^ u^^^ where u^^^ G • These indices 

j=0 j=0 

define the tensor products of spin matrices by S{u) = '^"ZqS^u). We also use the vector 
symplectic product introduced in Eq. (jl^ . When p = 2 we need the usual factor of —i if 
^(^■) = (1,1). 

The basic structure of the classes of indices defined by the mapping M is discussed in 
section EH and Appendix IX BI Specifically, class Ca of V2 (p") maps onto an n-dimensional 
subspace of V2„ (p). Each subspace is spanned by a set of n vectors Ga as defined in Eq. (jlH)) 
that depend explicitly on the parameters a = ( oq, ai. . . . , a„_i) in GF (p) which define a 
in GF (p") as a vector over GF{p). Since u o v = for any two vectors in Ca, it follows 
gr{ci) ° gs(tt) = for two generating vectors. 
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As in the case of n = 1 and n = 2, each non-zero vector in one of the Ca is mapped into 
a w 7^ G V2n{p) that can be written uniquely as 

n-l 

w = I]^igJ(")• 
J=o 

Assume p is odd. Following the paradigm established earlier, for a given density p define 

Xp(w) = Tr [p (r/-^go(a))'° " " " (^'""-^5g„-.(a))'""^] • (65) 
A discrete Wigner function for a density p on H^" is defined 

w 

where u o w is defined in Eq. (jl^ . 

It is not difficult to show that W^p(u) is real and Ylu^pi^) ~ '^^^ proof is simply a 
matter of keeping track of the various representations: 

p-1 



P 



1 

!2n 



a (6o,...,fen-l)7^(0,...,0) 



6o 



_l + lyTr 



Tr p 



-/ + 5^P,(uog,(«) + r,) 



This immediately confirms that Wp is real and shows that Wp (u) is the coefficient of the 
Hermitian matrix = (— / + J2a Poi{u o gj(a) + rj)) /p". 

For the normalization, summing over u is equivalent to summing over all of the vectors 
in each a summand: 



_pn + ^Tr {pi) 



as required. Again note that we inserted a factor of rj^'' into the gk{c() term to define a set of 
Wigner functions. This latitude of definition is exploited in the Appendix to give complete 
separabilty when p is an odd prime. Furthermore, with this special choice of phase factors, 
the analog of Eq. ()48|1 holds and the generalization of the argument for n = 2 gives 



n-l 



A{u) 



A{u 



0> 



j=0 
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whereu = 0^=0 

For p = 2 the same calculations apply provided factors of —i are included where required. 
However the methodology establishing separability fails for n > 2, and as far as we can 
determine the Wigner function as defined above does not respect separability. 



IX. DYNAMICS 

For completeness, we conclude with a discussion of Hamiltonian dynamics in the present 
context. Starting from the Heisenberg - von Neumann equation for a (i— dimensional system. 

'^^=i[p,H]=z{Hp-pH). (66) 

{h = 1) we obtain a closed form for the dynamics of either the Wigner function or the 
characteristic function when d = p, a prime. 

Let p denote an odd prime. The spin coefficients of a density p are defined by 

Su = tr {Sip) (67) 

so that 



1 

P = - 
p 



(68) 



Xp [mua) 



In defining the Wigner function, however, we emphasized the role of the characteristic func- 
tions xp {f^Ua) rather than the spin coefficients, and we also noted that one could add phase 
factors. For this discussion we use 

tr(^p{r,^-'-Su}f^) a^p 
tr [p [Sup]'^) a=p 

since the extra phase factors simplify the analysis. The same convention will be used for the 
Hamiltonian H. Of course, the spin function and characteristic function are simply related. 
Using Eq. (122), we obtain for y = {yo,yi) = mUa 

Sy = r]'-'''''Xp{-y), (69) 

and the phase factors allow us to avoid making a = p an exceptional case in . Thus (|US| 
becomes 

i^,,2-«™^^(_„)5,_. (70) 



" P 
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Since the spin matrices are orthogonal, it is easy to show that 



dt 

where 



^L{w,u)xp{-u) (71) 



L (w, u) = \h {u - w) {r,^-''^^- - r^^-^--) . (72) 

Equation (jH^I) enables one to convert ()72j) to describe the dynamics in terms of the spin 
coefficients rather than the characteristic functions. It is easy to check that L is a Hermitian 
operator indexed by V2 (p), so that (f7T|) can be solved in closed form. 

The evolution of the system can also be expressed in terms of the evolution of the Wigner 
functions. Using Eq. together with the results above, we avoid explicit use of the A [u] 
operators. Since the Wigner function is real, using Eq. (j21I)) we can write 

w w 

Then taking the time derivative, using ()72|) and then inverting (fTSj) gives 

dt 

where 



:Y,L{v,y)w,{y), (74) 



L {v, y) = - [v'^^yxH (2 {y - v)) - v'y^^xn (2 (v - y))] (75) 
p 

is Hermitian on V2 (p) . 

This representation works best when the density p evolves in the convex hull of the MUB 
projections. As an example when p = 3, let the Hamiltonian be 

H = u (^0,1 + 

and take p (0) to be -P{i,i) (0) = | [^0,0 + + ?7S'2,2] • Computing L and finding its spectral 
decomposition leads to the expression of p (t) in terms of MUB projections as 

p{t) = ^[(l + 2cos(c^t))P(i,i)(0) + (l + 2cos(cut + 27r/3))P(i,o)(l) 
+ (1 + 2cos(cjt + 47r/3))P(i,2)(2)]. 

In the special case of p = 2, the necessity of selectively introducing a factor of —i modifies 
the form of L. Any density p can be written as 



1 
2 

2 [Soft + ■so,i5'o,i + si^o'S'1,0 + ■si^iS'i^i] . 
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where so,i = 'itT'x, si,o = ttIz and = —iniy, and the m's are real with square sum less than 
or equal to 1. Defining the characteristic function as before, 



Xp {j, k) = Tr [p {aj,kSj,k)] , 



we find Xp i'^) equals the corresponding m and 



Sj,k = {-iy^'Xp {j,k). 



Working through the differential equation leads to a similar form: 



dXp i'v) 
dt 



i^L{v,y)xp (y) 



(76) 



u 



with L a Hermitian matrix given by 



(77) 



Thus the structure of L is similar to the p > 2 case but with powers of i rather than powers 
of 77 = —1. That difference makes the corresponding equation for the Wigner function more 
complicated, and we do not present it here. Our conclusion is that the discrete Wigner 

function is not particularly useful for studying the dynamics of a two-lcvcl system. 

A similar approach works for n systems, and we record the results for n = 2. For u =u^^^ © 
u^^^ = bogo (a) + bigi (a) set 



where for a 7^ we use go (a) = (1, j/oO) 0, yoi) and gi (a) = (0, yio, 1, t/n) . Recall that 
Z/oi — yw- One can then prove for all cases of a that 




Xp (-u) = X*p (u) 





(78) 



again for all a. 



34 



Recall the vector symplectic product 



uow = (n(°),7i«)o(wW,u;«)=^ 

k=0 

and for convenience set (u, u) = u^^^f'^ + u^q\^^\ Then 

u 

Using the analogous representation for the Hamiltonian, we have the analogue of ()7H): 

^^^ = .$:L(w,u)x,(-u), (80) 

u 

where 

L (w, u) = ^XH (-W + u) - 7^2-^—) (81) 

is Hermit ian on V4 (p). 

The derivation of the dynamics in terms of the Wigner functions follows almost word for 
word the pattern in the n = 1 case, since the computation of the Wigner function in terms 
of the characteristic function is symbolically identical. This time 

L (v, z) = 1 [v'^^-'xH (2v - 2z) - v'^°^XH (2z - 2v)] (82) 
is Hermitian on V4 (p) and 



V 



^J2l{^,z)W,. (83) 

z 

When p = 2 = n, we obtain a structurally similar result, although as before powers of i 
appear instead of powers of rj. Letting u^'^^ = (^u^q\u'^^^ and u = m^^^ © u^^\ 

Xp (u) = (-i)^"'"^ 



and 



The operator in the sum is Hermitian, and again the transformation to the Wigner function 
context does not seem to be particularly useful. 
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X. APPENDICES 
A. Finite fields 

Reference 

A finite field is a finite set of elements that contains an additive unit and a multi- 
plicative unit 1, that K is an Abelian group with respect to addition, K* = K — {0} forms 
an Abelian group under multiplication, and the usual associative and distributive laws hold. 
The simplest example of a finite field is the set of integers modulo a prime number p that 
is denoted by Zp = {0, 1, ■ ■ ■ ,p — 1}. If p is not prime there are elements that do not have 
inverses, for example the set Zl = {1,2,3} does not form a multiplicative group because 
2^ = mod4. 

It can be shown that if i^T is a finite field, then \K\, the number of elements in K, is 
p", the power of a prime. Fields with the same number of elements are isomorphic and 
are generically denoted as the Galois field GF {p^). A field containing p^ elements, n > 1, 
can be constructed using an irreducible polynomial / of degree n that has coefficients in 
GF{p) = Zp. Let 

/ (x) = Cn_iX""^ + . . . + CiX + Co 

be such a polynomial. Let A ^ GF{p) denote a symbolic root of / (x) = so that 

A" = - (c„_iA"-i + . . . + ciA + Co) . (84) 

It can be shown that each element in GF (p^) can be represented as 

a(A) = ^afeA'=. (85) 

fc=0 

Addition and multiplication proceed in the usual manner with the replacement of powers of 
A greater than n — 1 reduced by using Eq. (jSD). While the explicit representation depends 
on the choice of /, the theory guarantees different representations are isomorphic. 
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As an example, we saw in Section IVII that if n = p = 2, then / (x) = x'^ + x + 1 and 
GF{A) = {0, 1, A, A + 1}. For p an odd prime and = 2 we noted that elements of GF[p'^) 
could be written as j + A; A, where j and k are in GF{p) and f{x) = x^ ~ D with D a 
quadratic non-residue modp. 

In addition, there is a trace operation defined on GF (p") that is linear over GF (p) and 
that maps GF (p") to GF (p). Specifically, if Aq, • • • , A„_i denote the n distinct roots of /, 
then 

n-1 

tr (a (A)) = ^a{Xr). 

r=0 

The elements a in GF{p'^) can thus be viewed as a vector space over the field GF{p) with 
basis |A'^:0<A;<n}. A dual basis {gk (A) : < k < n} can be defined such that elements 
of GF{p"-) also can be written as a linear combinations of the g^s with coefficients in GF{p). 
The definition of a dual basis uses the trace operation with the requirement that 

tr [X^gk{X)]=Hj,k). 

This structure was described in the Appendix of and the complete theory is presented 
in MJ. 



B. Mutually unbiased bases for d=p". 

For the finite field GF{p^), as is explained in section IVTl we start with a vector space 
V2(p"). We need to map the vectors in V2(p") onto the space V2„(p) in order to write out 
the spin matrices corresponding to the set of MUB. A typical vector j3ua can be written as 

n-1 

(^u^ = J2 {x'^'\a,(3)e,+y'^^\a,(3)f,) . (86) 

j=0 

The x^^^ (a, (3) and y^-'-' (a, /?) are in GF (p) and {ej, : < j, < n} is a set of 2n linearly 
independent vectors over GF (p") . It is convenient to take them to be of the form ej = 
(1, 0) and fk = ^fc(A) (0, 1) so that 

triUoe,)=tr {X^ g^ (X)) = 6 {j, k) . 

The key point to defining a MUB is that for two non-zero vectors in V2 (p"), say 7iMa 
and 72tt/3, "jiUa ° 72^/3 = iff a = /5. Consequently, if in Eq. (|HB| we set Xr '' = x'^^^a, (3r) 
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and i/r''^ = y^^\a, (3r) for r = 1 and 2, we have 



= tr {piUa O P2Ua) 
/ n— 1 n— 1 



,j=0 A:=0 



n-1 



i=o 

n-1 

i=o 

Identifying the jth vector as the indices of the jth spin matrix in an n-fold tensor product, 
we have a necessary and 

sufficient condition for commutativity: 

(S)^—nS (j) (j) ®^_n S (k) (k) = <S)2—nS (k) (fe) ®?_n S (J) (j). 

Thus the set of vectors {'jUa, 7 G GF{p"')} corresponds to a commuting class 971^ of 
tensor products of spin matrices. The hnear mapping M : V2 (p") — > (p) defined by 



is one-to-one and onto. Using Eq. (j87j) this partitions the generahzed spin matrices into 
d + 1 commuting classes having only the identity in common and satisfying the condition 
for the existence of a set of c? + 1 mutually unbiased bases. In writing the M mapping we 
are using a different definition of the basis {cj, fj} that the one used in 2^. The definition 
in this paper lends itself more readily to a discusion of separability. 

C. Separability and the M mapping 

We provide some details about the mapping M : V2 (p") — > V2„ (p) . Let A denote a root 
of an nth order irreducible polynomial over GF{p). On V2{p^) recall the set of vectors 

{ej = A^(1,0), /, = g,iX)iO, 1), j = 0, 1, ■ ■ ■ , n - 1} 
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n-1 

where tr{fj o e^) = 6{j, k). Let a = '^j-^"' ^ GF{p"-) and using Eq. (jHUj) define 

i=o 

n-l 

Ua = {l,a) = eo + ^yf\a)fj 

j=0 

n-l 
fc=0 

Then for / = 1, • • • ,p — I 

AV = e, + 5^ yf\a)f) 

j=0 

where 

n-l 

l/f(«) = 5^tr(A^+'+'=)a,. (90) 

k=0 

Let us work out the details for the case p an odd prime and n = 2. We choose as our 
irreducible polynomial — D = mod p, where D is a quadratic non-residue. The symbolic 
roots of this equation are A and {p — 1)A. For example if p = 3 we may take D = 2. Then 
tr[/(A)] = /(A) + /(2A). It is not difficult to show ^o(A) = 2'^ and ^i(A) = {2D)-^\. Then 

y^^^ = 2ao, y?^ = y^^^ = 2Dai, ^ = 2Dao. 

We now can define the index generators of the MUB by 

G„ = |g.(a) = M(A^Mj = 0n(^)(a), ul^\a) = {6{j,r),yl^^), r = 0, 1, ■ ■ ■ , n - l| 

Gp. = |g.(p'^) = 0«(^)(p"), M(^nP") = (0,5(j,r)), r = 0,l,---,n-l|. (91) 

We should note that gr(p") is not M(A''(0, 1)) but rather M{gr{\){0, 1)). For the example 
of odd p and n = 2 we find for a = ao + aiX, qq, ai G GF{p), 

Ga = {go(a) = (l,2ao) © {0,2Da,), g^{a) = {0,2Da^) © (l,2Dao)} 
Gr,2 = {go(/) = (0,l)©(0,0),gi(p2) = (0,0)©(0,l)} (92) 

Each generator set is characterized by two independent four- vectors that determine a plane 
containing p"^ points. These planes intersect at only one point, the origin, and so the p"^ + 1 
sets determine p"^ — 1 distinct points and, including the origin, every point of V4{p). 
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We note from Eqs. (jH^ and that y'^^^ = which ensures the symplectic product 
is preserved by the mapping. Therefore, we have for the general case 

n— 1 n— 1 

e V,{pn ^ ui^\a) e V,^{p) ^ 5^„(„) = (g) Sl^, (93) 

j=0 j=0 

where ui''^ depends on a and the bj. With this notation, the mapping from the index space 
to the spin matrices is complete, 

{n— 1 n— 1 n—1 \ 

r=0 j=0 r=0 J 

For the case of an odd prime p and n = 2 this result is Eq. (jUj). 

The spin matrices can be further expanded with the help of Eqs. (jlOp and ()12p . the 
symmetry of the and a lot of algebra. First 

n—l n—1 n—1 

(8)n^!o)(.) = (8)^^...(")^*^^"'^ (94) 

j=0 r=0 j=0 
n-1 



r=0 

b) = (2-\{b, - l)yf{a) + 6, ^rZ/f («) ) • (95) 

\ r=0 

If 6j 7^ 0, define qj{a, b) = bjq'j^a, b) and we have 

C c6i -2-l6,(6,-l)Q;(a,, 

>Jbj,qj{a,b) — >Jl^q'Ja,b)'l 



If 6j = 0, we have 



iJO,gj{a,b) — iJo,! 



After some manipulation, we can then rewrite Eq. ()94|1 as 



n—1 n—1 



where the proper order of the tensor products is understood and where 

r jjLr bj=0 
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We now can incorporate the factor into the definition of x as is done in Eqs. ()51|) and ()65p. 
Again leaving the ordering of the tensor products understood, Eq. (jUBj) can be rewritten as 

n (')-"'='""'"'°'Ss,(o,)'' = (r"^s„,0» . (98) 

r bjT^O bj=0 

Therefore, we have shown that by introducing an appropriate phase factor that depends on 
r and a with each 5'g^(a) and by using rj^'^ ^ Sq^i in the definition of the one particle Wigner 
function, we can define a Wigner function for all n > 1 that respects complete separabilty 
for odd p. Note that the spin matrices appearing in the direct product are all in the standard 
form Su^ where c E Lp. 

For the example of odd prime p and n = 2 we have for ho and hi not equal to zero 

go(«,&) = ho2aQ + hilDai, qi{a,h) = hQ2Dai + hi2Dao 

e{a,h) = {ho + hi)2Dai. (99) 

As stated in Section IVII Dl the analysis for p = 2 requires special handling. For the case 
of a bipartite system, it was shown in IVII Dl that we could still prove a form of separability; 
however, for n > 2 we have been unable to make the method used here work. 



D. Symplectic structure of the MUB 

We have seen that Eq. (j^T|) determines the index sets for the MUB. If u,v G V2 {p^) , 
consider the transformations A : V2(p") V2(p") that leave the symplectic product u o v 
invariant. This is the set of 2 x 2 matrices wit^ 



which forms the symplectic group 5*^(2, K) 



1 entries in = GF(p'^) with unit determinant 

3, HQ 

We now want to study the mapping M defined in section IVII For simplicity we take n = 2 
so that the sets of generators of the MUB on Vi^p) are {go(a;), gi(a;)}. Introduce the 2x2 

/o-i\ f ^\ 

matrix o" = and the 4x4 matrix J = , then we can write the symplectic 

product in terms of an ordinary inner product, go(tt) ogi(a) = (go(a), Jgi(a)). Let A ^ A4 
where ^4 is a linear transformation on V4(p), such that A4go(Q;) o yl4gi(a) = go(«) o gi(a). 
Then a matrix representation of A^ must satisfy A4JA4 = J, where A/ is the transpose of 
A4. The set of linear transformations that satisfy this condition forms the symplectic group 
5*^(4, Zp). This is analogous to the canonical transformations for the continuous case. Under 
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such a transformation, the classes Ca determining the ONB of a given MUB are mapped into 
one another. In summary, the symplectic group 5*^(2, p") can be mapped onto a symplectic 
group Sp{2n, Zp) and the operators act on the bases in a MUB in such away as to leave 
the MUB invariant. For further discussion of the symplectic group in this context see Q]. 



E. Phase Space and Finite Geometry. 

The purpose of this section is to review the role played by the geometry of the phase space. 
In Section HVI we defined V2 {p) to be the phase space for the discrete Wigner function when 
n = 1, and lines in the vector space play an important role in relating the Wigner function 
to probability measurements. By analogy, for d = p"' a natural candidate for phase space for 
a d -level system is a two dimensional vector space with entries from an appropriate set of 
scalars which has d elements in it; that is, we consider V2 (p") = {{a, (3) : a, p e . 
However, in analogy with the continuous case for n each described on a Hilbert space Hp 
we use V2n{p) as the phase space. The M mapping takes lines in V2(p") to hj^erplanes in 

V2n{p)- 

If K denotes a finite field, the definition of a line in V2 {K) is the obvious one. A line L 
in V2 (K) is a set of points in V2{K) 

{{x, y) : -Xy + fix + = x,y e K} . 

We always omit the case in which A = /i = 0. It is important to note the line consists of 
these points and only these points. For example in V2(3), the sets Li = {(0, 0), (1, 1), (2, 2)}, 
L2 = {(0, 0), (1, 2), (2, 1)}, and L3 = {(0, 1), (1, 2), (2, 0)} are lines. Two lines intersect only 
if they have a point in common, otherwise they are parallel. The lines Li and L2 in the 
above example intersect at the origin, while Li and L3 are parallel. 

V2 (K) is also an example of an affine plane, a concept defined axiomatically in terms of 
a finite number of points, a finite number of lines, and the relationship that a point lies on 
a line. It can be shown that if a finite affine plane AP exists, then there is an m such that 
AP has exactly points, m? + m lines, each line contains m points and each point is on 
m + 1 lines. Two lines are said to be parallel if they have no point in common, and there are 
m + 1 sets of m parallel lines. (See for a summary of these results and references.) Since 
no affine plane is known for an m which is not a power of a prime, we are again restricted 
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to dimension p". 

The image under M of lines in V2 (p") , play a central role in the definition of a Wigner 
function, and we summarize a few of their properties. Using the generalization of Eq. (|31|1 
with GF{p^) replacing Zp we have: 

L{a, 7) = {xUa + JUd : x G GF^p"^)} 
^(p",7) = {yud + iuo:yeGF{p^)} 

with a,7 G GF(j/'). Recall that = and Up^ = (0,1). The vectors Ua and Ud 

multiplying the variables x or y were introduced earlier in Eq. (j39|) as a convenience. They 
now are playing the role of "slopes" in an indexing of lines in V2 (p^), a much more general 
setting. For each slope, as 7 varies over GF{p'^) we get a set of parallel lines that contains 
each point in V2(p") once. 

1. Each line contains elements, and there are p^" +p" distinct lines. 

2. The lines through the origin, L(a,0) where a & Ipn = GF{p"') U {p"}, only intersect at 
the origin. Furthermore, 

U (L(a,0)-{(0,0)}) = \/2(p")-{(0,0)}. 

3. Each set of parallel lines partitions V2{p^): 

V,{p-)= U L(a,7) 

7eGF(p") 

for each a G Jpn. 

The relevance to this paper of the affine plane is that it can be shown that for certain 
values of d, such as c? = 6, there is no corresponding affine plane, and for other values of 
d, such as (i = 12, the existence of a corresponding affine plane is an open question. (See 
standard texts in combinatorics for more details or [8] for references.) We have already 
noted that if K denotes a finite field, then V2 (K) is an example of an affine plane, so we 
are working in the most general context with the necessary structure. 

The last tool we need is the symplectic product of vectors in V2 (K) over the finite field 
K. Specifically, recall that 

(/il, Z/i) O {fi2, 1^2) = ^lfi2 - fil^2, (100) 



43 



where the algebra is in the field K. As an example, o = if and only M a = j3. For 
each a ^ d 'm. Id, (A, yu) is on the line L{a, 7) where 7 = (A, /i) o Ua- 

Finally, from each we generate n linearly independent vectors that are mapped into 
an n-dimensional hyperplane in V2n{p) using Eqs. (j89j) and (j9Up . 



F. Examples of the geometry 

1. One qubit 

Let K = GF (2), the Galois field consisting of the integers mod 2. The six lines of V2 (2) 
fall into three classes containing two parallel lines: 

{L(0, 0) = {(0, 0) , (1, 0)} , L(0, 1) = {(0, 1) , (1, 1)}} 

0) = {(0, 0) , (1, 1)} , 1) = {(0, 1) , (1, 0)}} 
{L(2, 0) = {(0, 0) , (0, 1)} , L(2, 1) = {(1, 0) , (1, 1)}} 

2. Two qubits 

The elements of A' = GF (2^) can be represented as 

{0,1,A,A2 = A + 1}, 

where is the additive identity, 1 is the multiplicative identity, and 1 + 1 = A + A = 0. The 
other relations follow in the obvious way, such as A (A + 1) = A^ + A = A + 1 + A = 1. The 20 
lines of V2 {GF (4)) fall into five classes of four parallel lines each. The class of vertical lines 
is generated by L(4, 0) = {(0, 0) , (0, 1) , (0, A) , (0, A + 1)}, and shifts of L(4, 0) by 7 (1, 0) . 
The other four classes are generated by L{a, 0) = {Pua '■ P G GF (4)} and shifts by 7 (0, _1)_, 
where L(0, 7) corresponds to a horizontal line. Graphs of lines in V2 (2^) appear in both 
and ba. 



G. Positivity relation 

We include this brief discussion in order to illustrate the difficulty in determining whether 
a given phase space function corresponds to a positive operator. The method given here is 
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closely related to the proof given in [l6]. Let {cjk} be an arbitrary set of complex coefficients 
and define the matrix B = 'Ylij k^i^Sj^k- Then p > if and only if tr{pBB'^) > for all B. 
Writing out the sum and using the properties of the spin matrices gives 

tripBB^) = J2 c,kc:MpS,-s,k-t)v'^'-''^- 

j,k,s,t 

Now we can express the trace in terms of the characteristic function 



xiyup) if x = 
x{xUa) where a = x~^y, x 7^ 



Therefore, we have, a not very illuminating, necessary and sufficient condition for x to arise 
from a positive matrix. The necessary and sufficient condition for x to correspond to a 
density matrix also requires that x(0) = trp = 1. 



[1] S. Bandyopadhyay, P. O. Boykin, V. Roychowdhury, F. Vatan, Algorithmica 34, 512 (2002), 

Iquant-ph /OlOS]^ (Sept. 2001). 
[2] P. Bianucci, C. Miquel, P. Paz, and M. Saraceno, | quant -ph/0106091 , (Jan. 2001). 
[3] J.S. Bell, Physics 1,195 (1964). Reprinted in J. S. Bell, Speakable and unspeakable in quantum 

mechanics, Cambridge University Press, Cambridge, 1987. 
[4] A. R. Calderbank, P. J Cameron, W. M. Kantor, J. J. Seidel, Proc. London Math. Soc.3, 

436-480 (1997). 
[5] D. Fivel, Phys. Rev. Lett. 74, 835 (1995). 

[6] G. B. Folland, Harmonic Analysis in Phase Space, Princeton University Press, Princeton, NJ, 
1989, W. P. Schleich, Quantum Optics in Phase Space, (Wiley- VCH, Berlin, 2001). 



[7] E. F. Galvao, |quaiit-ph/0405 070 , (May 2004). 

[8] K. S. Gibbons, M. J. Hoffman,W. K. Wootters, quant-ph/040115 v3 (2004). 

[9] M. Hillery, R. F. O'Connell, M.O. Scully and E. P. Wigner, Phys. Rep. 106, 121 (1984). 

[10] L. K. Hua, Introduction to Number Theory, (Springer, New York, 1982), p. 364. 

[11] I. D. Ivanovic, J. Phys. A 14, 3241 (1981). 

[12] M. Koniorczyk, V. Buzek, and J. Hanszky, Phys. Rev. A 64, 034301 (2001). 

[13] U. Leonhardt, Phys. Rev. A 53, 2998 (1996). 



45 



[14] R. Lidl and H. Niederreiter, Finite Fields, Encyclopedia of Mathematics and its Applications 
Volume 20, Addison- Wesley, Reading MA, (1983). 



[is: 

[16 

[17: 

[18 
[19 
[20 

[21 
[22: 

[23 
[24 
[25' 
[26 
[27 
[28 
[29 



Longenacker and Roettler, quant-ph/0309120 (Sept. 2003). 

F. J. Narcowich and R. F. O'Connell, Phys. Rev. A 34, 1 (1986). 

Juan Pablo Paz, Augusto Jose Roncaglia, Marcos Saraceno, quant-ph /04 1 1 1 7{ (Oct. 2004). 
A. Peres, Phys. Rev. Lett. 77, 1413 (1996). 

A. O. Pittenger and M. H. Rubin, Phys. Rev. A 62, 032313, (2000). 

A. O. Pittenger and M. H. Rubin, Linear Alg. Appl. 390, 255 (2004), |quant-ph/0308142l 
(Aug. 2003). 

J. Schwinger, Proc. Nat. Acad. Sci. 46, 570 (1960). 

A. Terras, Fourier Analysis on Finite Groups and Applications, Cambridge University Press, 
Cambridge, 1999. 

J. A. Vaccaro and D. T. Pegg, Phys. Rev. A41, 5156 (1990). 
A. Vourdas, Rep. Prog. Phys. 67, 267 (2004). 

H. Weyl, The theory of groups and quantum mechanics, (Dover Publications, New York, 1950). 

E. P. Wigner, Phys. Rev. 40, 749 (1932). 

W. K. Wootters, Found. Phys. 14, 391 (1986). 

W. K. Wootters, Ann. Phys. 176, 1 (1987). 

W. K. Wootters and B. D. Fields, Annals of Physics 191, 363-381 (1989). 



46 



